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Abstract 

An important problem in the theory of compressible gas flows is to understand the 
singular behavior of vacuum states. The main difficulty lies in the fact that the system 
becomes degenerate at the vacuum boundary, where the characteristics coincide and 
have unbounded derivative. In this paper, we overcome this difficulty by presenting a 
new formulation and new energy spaces. We establish the local in time well-posedness 
of one-dimensional compressible Euler equations for isentropic flows with the physical 
vacuum singularity in some spaces adapted to the singularity. 

1 Introduction 

One-dimensional compressible Euler equations for the isentropic flow with damping in Eule- 
rian coordinates read as follows: 

Pt + {pu) x = 0, 
pu t + puu x +p x = -pu, 

with initial data p(0,x) = po(x) > and u(0, x) = uq(x) prescribed. Here p, u, and p 
denote respectively the density, velocity, and pressure. We consider the polytropic gases: 
the equation of state is given by p = Ap' 1 , where A is an entropy constant and 7 > 1 is 
the adiabatic gas exponent. When the initial density function contains vacuum, the vacuum 
boundary T is denned as 

T = cl{(t, x) : p(t, x) > 0} n cl{(t, x) : p(t, x) = 0} 

where cl denotes the closure. A vacuum boundary is called physical if 

dc 2 

0<|^|<oo (2) 



in a small neighborhood of the boundary, where c = y j^p(p) is the sound speed. This phys- 
ical vacuum behavior can be realized by some self-similar solutions and stationary solutions 
for different physical systems such as Euler equations with damping, Navier-Stokes equations 
or Euler- Poisson equations for gaseous stars. For more details and the physical background 
regarding this concept of the physical vacuum boundary, we refer to [10\ [T9] . 

Despite its physical importance, even the local existence theory of smooth solutions fea- 
turing the physical vacuum boundary has not been completed yet. This is because the hyper- 
bolic system becomes degenerate at the vacuum boundary and in particular if the physical 
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vacuum boundary condition ([2]) is assumed, the classical theory of hyperbolic systems can 
not be applied [10] : the characteristic speeds of Euler equations are u ± c, thus they become 
singular with infinite spatial derivatives at the vacuum boundary and this singularity creates 
a severe analytical difficulty To our knowledge, there has been no satisfactory theory to treat 
this kind of singularity. The purpose of this article is to investigate the local existence and 
uniqueness theory of regular solutions (in a sense that will be made precise later and which 
is adapted to the singularity of the problem) to compressible Euler equations featuring this 
physical vacuum boundary. 

Before we formulate our problem, we briefly review some existence theories of compressible 
flows with vacuum states from various aspects. We will not attempt to address exhaustive 
references in this paper. First, in the absence of vacuum, namely if the density is bounded 
below from zero, then one can use the theory of symmetric hyperbolic systems; for instance, 
see [12] . In particular, the author in [T7] gave a sufficient condition for non- global existence 
when the density is bounded below from zero. When the data is compactly supported, there 
are two ways of looking at the problem. The first consists in solving the Euler equations in 
the whole space and requiring that the equations hold for all x and t £ (0, T). The second 
way is to require the Euler equations to hold on the set {(t,x) : p(t,x) > 0} and write an 
equation for the vacuum boundary T which is a free boundary. 

Of course in the first way, there is no need of knowing exactly the position of the vacuum 
boundary. The authors in [15] wrote the system in a symmetric hyperbolic form which allows 
the density to vanish. The system they get is not equivalent to the Euler equations when the 
density vanishes. This was also used for the Euler-Poisson system. As noted by the authors 

7 — 1 

\13\ 114], the requirement that p~^~ is continuously differentiable excludes many interesting 

solutions such as the stationary solutions of the Euler-Poisson which have a behavior of the 

i 

type p ~ \x\ t- 1 at the vacuum boundary. This formulation was also used in [2] to prove the 
local existence of regular solutions in the sense that p 1 ^ - ,u G C([0, T); H m (M. d )) for some 
m > 1 + d/2 and d is the space dimension (see also [16] . for some global existence result of 
classical solutions under some special conditions on the initial data, by extracting a dispersive 
effect after some invariant transformation, and [3]). 

For the second way and when the singularity is mild, some existence results of smooth 
solutions are available, based on the adaptation of the theory of symmetric hyperbolic system. 
In [9] , the local in time solutions to Euler equations with damping (p} were constructed when 
c a , < a < 1 is smooth across T, by using the energy method and characteristic method. 
They also prove that C 1 solutions cross F can not be global. However, with or without 
damping, the methods developed therein are not applicable to the local well-posedness theory 
of the physical vacuum boundary. We only mention a result in [18] for perturbation of a 
planar wave. For other interesting aspects of vacuum states and related problems, we refer 

to uning. 

As the above results indicate, there is an interesting distinction between flows with damp- 
ing and without damping, when the long time behavior is considered. Indeed, without damp- 
ing, it was shown in [8] that the shock waves vanish at the vacuum and the singular behavior 
is similar to the behavior of the centered rarefaction waves corresponding to the case when c 
is regular [10] . On the other hand, with damping, it was conjectured in [7j that time asymp- 
totically, Euler equations with damping ([1]) should behave like the porus media equation, 
where the canonical boundary is characterized by the physical vacuum boundary condition 
([2]). This conjecture was established in [3] in the entropy solution framework where the 
method of compensated compactness yields a global weak solution in L°° . We point out 
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that the difficulty coming from the resonance due to vacuum in there is very different from 
the difficulty that we are facing, since we want to have some regularity so that the vacuum 
boundary is well-defined and the evolution of the vacuum boundary can be realized. 

In order to understand the physical vacuum boundary behavior, the study of regular 
solutions is very important and fundamental; the evolution of the vacuum boundary should be 
considered as the free boundary/interface generated by vacuum. In the presence of viscosity, 
there are some existence theories available with the physical vacuum boundary: The vacuum 
interface behavior as well as the regularity to one-dimensional Navier-Stokes free boundary 
problems were investigated in And the local in time well-posedness of Navier-Stokes- 
Poisson equations in three dimensions with radial symmetry featuring the physical vacuum 
boundary was established in [5]. On the other hand, the free surface boundary problem 
was studied in [6] in the motion of a compressible liquid with vacuum by using Nash-Moser 
iteration; the physical boundary was treated in a sense that the pressure vanishes on the 
boundary and the pressure gradient is bounded away from zero, but the density has to be 
bounded away from vacuum, and thus the analysis is not applicable for the motion of gas 
with vacuum. 

In the next section, we formulate the problem and we state the main result: The vacuum 
free boundary problem is studied in Lagrangian coordinates so that the free boundary be- 
comes fixed. By change of variables, the equations can be written as the first order system 
with non-degenerate propagation speeds that have different behaviors inside the domain and 
on the vacuum boundary. In order to cope with these nonlinear coefficients, which give rise 
to the main analytical difficulty, the new operators V, V* are introduced. Our theorem is 
stated in V, V* framework. 

2 Formulation and Main result 

We study the initial boundary value problem to one-dimensional Euler equations with or 
without damping for isentropic flows ([T]). First, we impose the fixed boundary condition on 
one boundary x = b : u(t,b) = 0. The class of the initial data po, uq of our interest is 
characterized as follows: for a < x < b, where — oo < a < 6 < oo 



The condition (i) implies that the initial vacuum is physical, (ii) means that x = a is the 
only vacuum, (Hi) represents the finite total mass of gas, (iv) is the compatibility condition 
with the boundary condition at x = b. We seek p(t,x), u(t,x), and a(t) for t G [0, T], T > 
and x 6 [a(t), b], so that for such t and x, 



(i) p {a) = , < — p^ \ x=a < oo ; 
(ii) po(x) > for a < x < b ; 



(Hi) / po(x)dx < oo ; (iv) uo(b) = . 




p(t, x) and u(t, x) satisfy ([T]) ; 
p(t,a(t)) = 0; u(t,b) = 0; 
d 

< 7^/ j7 " 1 U'=a(t) <°°- 
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For regular solutions, the vacuum boundary a(t) is the particle path through x = a. In one- 
dimensional gas dynamics, there is a natural Lagrangian coordinates transformation where 
all the particle paths are straight lines: 

y = I p(t,z)dz, a(t) < x < b . 

Ja(t) 

Note that < y < M, where M is the total mass of the gas. Under this transformation, 
the vacuum free boundary x = a(t) corresponds to y = 0, and x = b to y = M; thus both 
boundaries are fixed in (t,y). By this change of variables, 

d t = d t - ud y , d x = pd y 

the system ([1]) takes the following form in Lagrangian coordinates (t,y): for t > and 
< y < M, 



(3) 

Ut+Py = ~U 



p t + p U y = 



where p = Ap 1 with 7 > 1. The boundary conditions are given by p(t, 0) = and u(t, M) = 
0. The physical singularity ([2]) in Eulerian coordinates corresponds to < \p y \ < 00 in 
Lagrangian coordinates and thus the physical vacuum boundary condition at y = can be 
realized as 

p ryj yi for y ~ . 



Euler equations ([3|) can be rewritten as a symmetric hyperbolic system 

(f> t + pUy = , ^ 
U t + H4>y = —U , 

in the variables 



q> = p 2 and u = \f A^p 2 . 

7-1 

Note that the propagation speed p becomes degenerate and the degeneracy for the physical 

7+1 

singularity is given by p ~ y 2 ~t . In order to get around this difficulty, we introduce the 
following change of variables: 

27 7-1 7+1 
£ = -y 2 ~> such that d y = y 2 ~t . 

We normalize A and M appropriately such that the equations ([!]) take the form: 

i.iil 

<k + (j)~<- 1 H = 0, 

(h 7+1 
u t + (^r) 7 - 1 ^ = ~u, 

for t > and < £ < 1. The physical singularity condition < \p y \ < 00 is written as 
< \<p^\ < 00. Thus we expect <p to be more or less £ for short time near 0. Since the 
damping is not important for the local theory, for simplicity, we consider the pure Euler 
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equations. Letting 

*-t T = il±i. 

7 2 7 - 1 ' 

the Euler equations read in (t, £) as follows: for i > and < £ < 1, 

^ + (|) 2 S = o, 

n t + (|) 2 V ? = 0. 

The range of is ^ < A; < 00, since 7 > 1. When 7 — ► 1, fc — ► 00 and when 7 = 3, we 
get k = 1. Note that the propagation speed is now non-degenerate. However, its behavior is 
quite different in the interior and on the boundary, since lim^o ^ = (0), but 4> < ^ < c<p 

if £ > ^ > 0. This makes it hard to apply any standard energy method to construct solutions 
in the current formulation. 

We will propose a new formulation to (|5|) such that some energy estimates can be closed 
in the appropriate energy space. As a preparation, we first define the operators V and V* 
associated to (0 as follows: 



1 6 2k 6 2k 1 



for /, g £ L 2 -, where we have denoted L^[0, 1] by Lg. We can think of V and V* as modified 
first order spatial derivatives. We also incorporate the boundary condition at £ = 1 in the 
domain of of V and V*, namely V and V"*, are given as follows: 

V(V) = {/ € L\ : e 2$ 

= {g G L| : G L|, 5 (£ = 1) = 0.} 

We also introduce the higher order operators (V) 1 and (V*) 1 : for / G T>(V) and g G P(V*), 

po , (/) = /tvvyv) i» = 2i 

|V(V*Vp(/) ifi = 2j + l 

and the associated function spaces X k ' s and y fe ' s for s a given nonnegative integer: 

X k > s = {/ G L\ : G Lf, < i < s} 

Y k > s = {geL 2 : {V*)\g) G L 2 , < i < s} 

equipped with the following norms 



(9) 



i=0 " i=0 

In order to emphasize the dependence of k, equivalently 7, we keep k in the above definitions. 
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In terms of V and V*, the Euler equations (|5|) can be rewritten as follows: 



i t ( io ) 

with the boundary conditions 

<f>(t,0) = and u(t, 1) = 0. (11) 

In this new V, formulation, the system is akin to the symmetric hyperbolic system with 
respect to V,V*. In particular, the zeroth energy estimates assert that this V, V* formula- 
tion retains the energy conservation property, which is equivalent to the physical energy in 
Eulerian coordinates: It is well known that the energy of Euler equations without damping 
is conserved for regular solutions: 



-{/' - 

dt Ja(t) 2 



pu 2 H -dx} = 0. 

(t) 2 7-1 



and in turn, since dy = pdx, it is written as, in y variable, 

d f M 1 A 
dt Jo 2 7-1 

It is routine to check that this is equivalent to 

1 d . f 1 2+1 o 7 — 1 +1+1 ,9 , -, 

5*U (-- 2 + 7 —^-'^}-^ 

which is exactly the zeroth energy estimates of (|10[) with respect to V, V*: 

-{^ 1 ^ T leV| 2 + l^n| 2 ^} = 0, (12) 



ld_ 
2~dt 



This verifies that this V, V* formulation does not destroy the underlying structure of the 
Euler equations. Indeed, it also captures the precise structure of the singularity caused 
by the physical vacuum boundary, and furthermore, it enables us to perform V, V* energy 
estimates yielding the a priori estimates and the well-posedness of the system. 

In order to state the main result of this article precisely, we now define the energy func- 
tional £ k ' s (4>, u) by 



^••(0.«) = 5iTTTll^lfe + llf*«lli? 



1 

2* + l- — i ' - (13) 

+ {2k+l)2 \\v(em 2 Y k ,^ + \\v*(eu)\\ 2 xk ,s-i ■ 

To close the energy estimates, s will be chosen as s = \k~\ +3, where \k] is a ceiling function, 
namely \k] = min{n G Z : k < n}. We are now ready to state the main results. 

Theorem 2.1. Fix k, where \ < k < oo. Suppose initial data 4>q and uq satisfy the following 
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conditions: 

(i) £ fc >W+3( 0O)Uo ) < _L < ^ < c for some C > 1 

There exist a time T > onZy depending on £ k ^ k ^ +3 (4>o,uo) and Cq, and a unique solution 
((/), u) to the reformulated Euler equations (|10|) wrei/i i/ie boundary conditions (jlip on t/ie time 
interval [0, T] satisfying 

and moreover, the vacuum boundary behavior of (j) is preserved on that time interval: 

1 6 



Remark 2.2. The evolution of the vacuum boundary x = a(t) is given by a(t) = u(t,£ = 0). 
By Theorem \2.1l one can derive that u^(t,^) is bounded and continuous in (t,x) € [0, T] x 
[0,1], and since u(t,0) = f? u^(t, £)<i£, we deduce that the vacuum interface is well-defined 
and within short time t <T, the vacuum boundary at time t stays close to the initial position 
with \a(t) — a\ < CT for some constant C depending on the initial energy in Theorem \2.1\ 

Remark 2.3. Different constants in the energy functional f|13|) are due to the nonlinearity 
of ([5]) or (fTOl) . The structure of the equations are to be systematic after applying V,V* to 
(|10p as in (|40p and thereafter. Since we are dealing with the local existence theory, one may 
work with the energy functional £ k,s = | |^ fc , s 4- | \£ k u\ \yk, a , which is equivalent to (TL3j) . 

Our work is a fundamental step in understanding the long time behavior of regular solu- 
tions and the vacuum boundary of Euler equations with the physical singularity rigorously 
With or without damping case, it would be interesting to study whether our solution exists 
globally in time. In particular, for the damping case, the physical singularity is expected 
to be canonical as in the porus media equation, it would be also interesting to investigate 
the asymptotic relationship between our solution and regular solutions to the porus media 
equation. 

Parallel to the recent progress in free surface boundary problems with geometry involved, 
we expect that our result can be generalized to multidimensional case, since the difficulty of 
the physical singularity lies in how the solution behaves with respect to the normal direction 
to the boundary. We leave this for future study. The physical vacuum boundary also naturally 
appears in Euler- Poisson equations for gaseous stars. It would be very interesting to study 
the behavior of solutions and the vacuum boundary under the influence of the gravitation. 

The method of the proof is based on a careful study of V, V* operators and V, V* energy 
estimates. The first key ingredient is to establish the relationship between the multiplication 
with |, a common operation embedded in the equations ([5]) or (fTUj) . and V, V* by using the 
underlying functional analytic properties of V, V* which can be obtained from a thorough 
speculation on the behavior of V, V* near the vacuum boundary £ = 0. Another essential 
idea is to find the right form of spatial derivatives of df(f), which is critical to cope with strong 
nonlinearity, in particular of the second and third order equations and to close the energy 
estimates in the end. For large k, this can be done by introducing the representation formula 
of (V*) l (!; k u). In the similar vein, due to the strong nonlinearity, the approximate scheme 
starts from the third order equation and lower order terms including (f> and u are recovered 
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by integrals and boundary conditions. Lastly, we point out that it is not trivial to build the 
well-posedness of linear approximate systems and the duality argument is employed for that 
as in pp. 

The rest of the paper is organized as follows: In Section El we study the operators V 
and V* built in the reformulated Euler equations ([5]) or (|10|) . In Section HI we establish the 
a priori estimates in V, V* formulation. Based on the a priori estimates, in Section [5j we 
implement the approximate scheme and prove that each approximate system is well-posed. 
In Section^ we finish the proof of Theorem 12. 11 In Section[7l the duality argument is proved. 



3 Preliminaries 

Throughout this section, we assume that <fi is a given nonnegative, smooth function of t and 
£, and moreover, j and d^cj) are bounded from below and above near £ = 0. 

3.1 Basic properties of V and V* 

In this subsection, we study the operators V and V*. Let us denote C^°((0, 1)) (respectively 
Cq°((0, 1])) the set of C°° functions with compact support in (0, 1) (respectively (0, 1]). 



Lemma 3.1. 



-x 



C °°((0,1]) =V(V) = X k ' 1 ; 



-Y k ' 1 



Cg°((0,l))' = V(V*)=Y k > 1 ; 
Proof. Take / G X k >\ Hence, 

:m = J q ^(^/)| 2 + |/| 2 ^<oo. 



(14) 



We make the change of variable y = £ 2fc+1 and F(y) = ^rf. Hence 



2 



rl j f 4fc 1 

1=j / o (2* + l)ft^ + — i s - 5 -|FNv<oc 



Since, 2k > 1, we deduce that ^ttt > 1 and hence, necessary F(0) = 0. Applying Hardy 



2fe+l 

inequality to F, we deduce that 



l F 2 
—dy<C 

o y 



f\dyF\ 2 - 
Jo 



Hence, going back to the original variables, we deduce that for / G X k ' 1 } we have 

1 ^di<C\\i\\f T \\V{f)\\l r (15) 



Now, consider a cut-off function x G C°°(M) given by x(0 = f° r £ — 1/2 and x(£) = 1 
for £ > 1. We also define x n (C) = For / G A"*' 1 , we define /„(£) = Xn(0/(0- Hence, 
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f n G X k ' 1 and it is clear that / n goes to / in I? . Moreover, 



V(/ - /„) = (1 - XnM/) " nx'K)^/ • 



The first term on the right hand side goes to zero in L| when n goes to infinity. For the 
second term, we use the fact that 



L 



1 ±2k i p\ £% 



which goes to zero when n goes to infinity in view of Q15]). Hence, we deduce that / n goes 
to / in Now, it is clear that /„ can be approximated in X k > 1 by functions which are in 

Cq°((0, 1]). Indeed, one can just convolute /„ by some mollifier. Hence, the first equality of 
CHD holds. 

To prove a similar result for V*, we take g G Y k ' 1 and fix some c € (0, 1]. For £ G (0, c] , 
we have 

, r»f9tf)^t, , 5(c), 



< e ^ + ^ c 



where e = C J Q C l^(^^)| 2 ^ 7 - By choosing c small enough, we can make e small. 

Hence, we deduce that for all e > 0, there exists a constant C € (g) = |^r| such that for all 
£ G (0, 1], we have 

We denote g n = x n <7> hence it is clear that g n goes to g in L|. Moreover, 
The second term on the right hand side satisfies 



V*(g - g n ) = (1 - Xn)^*(<?) " n X '(nO^:g. 



JO 4 JO ? 2n -■»-». 

Hence, since e > is arbitrary, it goes to zero when n goes to infinity. Therefore, we deduce 
that g n goes to g in Y ,1 . Now, it is clear that g n can be approximated by functions which 
are in Cq°((0, 1]) and the second equality of (fH|) follows. □ 

Lemma 3.2. (%) For f G JC*- 1 , 5 G F^' 1 
(2) ker V = {0} and ker = {0} 
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(3) Vt andVj* are commutators ofV,V* and dt: 

d t V(f) = V(d t f) + V t (f), d t V*(g) = V*(d t g) + V t *(g) 

where 

In addition, 

V t *(g) = -2k^V*(g) and VV t *(g) = V t V*(g). 

Proof. The proof of this lemma directly follows from the density property proved in the 
previous lemma. It will be used frequently during the energy estimates. □ 

The following lemma displays a key ingredient of the main estimates. Dividing by £ is 
a common operation embedded in the equations ([5]) and (llOp , The lemma claims that the 
operation, which acts like derivatives when £ approach 0, is completely controlled by modified 
derivatives V and V* . 



Lemma 3.3. (1) If f G X k ' 1 and g G Y k ' 1 , then £ and | are bounded in and we obt 
the following inequalities: 



am 



< C\\^\\f ? \\Vf\\ LV < C7{||||||J.||^||^} (16) 

(2) More generally, if f G L| satisfies ^m-i G £| / or some nonnegative real number m, 
m + k > | , i/iera 

II ^ II <r- /^ll^ l|2fc ii II 

— \\t2 (_/ — roo — 7 r2 . 

1 tm 11 ^ — a) £ f m— ^ £ 

(3) Also, if g satisfies G L| and j^ff G L| /or some nonnegative real number m, 
m < k + i , i/ten 

II 9 II <r ril^l|2fc ||Z!#J| i ii 9 ii 
Here, m is not necessarily an integer. In practice, m will be chosen as ^ < m < k 



Proof. The point (1) was already proved for / in the previous lemma by Hardy inequality 
(see (fT5j) ). For the second inequality, consider first g G C^°((0, 1)), hence 



2 



and g(l) = 0. Since k ^ \ and 

l\r* n ^ 9 ac ^ ii £ n2fe in/*„n ii^ii 

we obtain the desired result. The case where g G Y k)1 follows by density. 

Now, we concentrate on (2). For <fi satisfying the same type of bounds as <j), we define 
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Va(f) = ^*<%(^r/)- For / as in (2), we use that 

V(f) _ 1 ^ , ^ f_ S = y f f v 

where cf> is given by ( ^ 2 ( fc + m ^ 1 ) = <^ 2fc £2(m— l) _ gi nc6) ^_|_ m _ j > 1 we can apply the estimate 
([15]) for V replaced by Vk+m~i and / replaced by . This gives the desired bound. 
For (3), we write 

V*(g) <t> 2 \~ r , 9 , 



-V, 



with </> = £. Since, m — 1 — k < — |, V m -\-k satisfies the same estimates as V*, in particular 
the second estimate of (fT6j) holds for V m -i-k, hence (fTT|) holds. □ 



Remark 3.4. We note that the boundary conditions on f and g at £ = in Lemma \3.3\ are 

imbedded in X k,s and Y k)1 , namely L| boundedness ofVf and V*g forces f and g to vanish 
at £ = 0. Actually, it forces them to be less than £ 1//2 . 

As a direct result of Lemma 13.3} we obtain the following L| estimates for ^4- and ^ m . 

Corollary 3.5. For given nonnegative real number m, m < k + |, if f 6 X k '^ m ^ , then there 
exists C\ only depending on so that 

W-J^Wl* < Cl||/llx*.M • 

Also, for given nonnegative real number m, m < k + \, if g G Y k '^ m \ then there exists C2 
only depending on so that 

ll 9 11 ^ n 11 11 

Proof. We apply (2) and (3) of Lemma [3.3l alternatinglv until negative powers of £ disappear. 
Note that the conditions on m come from the one in (3) of Lemma 13.31 □ 

We next prove the Sobolev imbedding inequalities of V, V* version, which will be useful 
tools to control nonlinear terms. 

Lemma 3.6. If f E X k ^ k ^ +1 and g G Y k '' k ' +1 , then there exist constants C3 and C4 only 
depending on ||^||z,|° so that 

f 9 
Ht^IU^ < C s \\f\\ X k,[k^+i and ||7jf||i«= < C 4 | \g\ \ Y k,m+i ■ 

Proof. We start with g part. By the definition of V and V* , one finds that 

8l 9_ ] = _e^Y*g 

?L^fcJ i 2 k £k 
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Thus, by Sobolev embedding theorem in one dimension, it suffices to show that 



9 V*g 



6L2. 



£ k ' £fc fc £ 

This follows from Corollary 13.51 Hence, we conclude that LT bound of S% is controlled by 



||g , ||yfe,rfcl+i. For / part, we show that |||2F^-||i°° is bounded by ||/||x fc .rfcl+i • Note that 

U tl£2k£kt £k ^ £2k £k+l - 

Thus, by applying Corollary 13.5} we obtain the desired conclusion. □ 

More generally, we obtain the following: 

Lemma 3.7. Let < j < k — ^ be a given nonnegative number. If f E X k, ^~^ +1 and 
g £ Y k '^^ +1 , then there exist constants C§ and Cq only depending on ||4||l°° so that 



IItjIIli < C 5 ||/|| X fc,rji+i and < Celbllyfc.m+i 

Proof. We only treat |j. Note that 



9^} = -^^ + ^ -3 



£J J 02fc £j V ■*> £j+\ 

Hence, by the Sobolev embedding theorem, it suffices to show that |y, and |^rr are in 
L|. This follows from Corollary 13.51 Note that j has to be less than k — ^. □ 

Next we present the product rule for the operators V, V* . 

Lemma 3.8. Let f,g€ DVn'DV* be given. Let h be a given smooth function. The following 
identities hold: 



g _ , 2k V{e<t>)9 



V ° = - V *S + 2k ^ 9 ^ = - V *° + 2k + l ^ 
V(fh) = V(f)h + f^d^h, V*(gh) = V*(g)h - g^d^h 



6 2k 6 Ak f k 6 2k 1 dp 6 

• ^[fg] = V(f)g + ^fd^g] = V{f)g - fV*(g) + 2k^(- - -f)fg 

The lemma tells that when V or V* act on a function h, depending on that function, they 
can yield | or V ^ k ^ | besides Vh or V*h. 



12 



3.2 Homogeneous operators V, V 

Next we introduce homogeneous linear operators V and V* of V and V* as follows: 



V(f) = ^d^ k f}, V* (g) = and g(£ = l)=0. (18) 



These homogeneous operators are special cases of V and V* for which (j) is simply taken as 
£. Function spaces X ' and Y ' for s a given nonnegative integer are given as follows: 

X KS = {/ G L\ : G L|, < < < s} 

Y k ' S = {g G L\ : (F*)^) G L|, < » < s} 

where (V)* and {V*) % are defined as in ([7]) and (JSJ). These spaces are equipped with the 
following norms 

||/||^^||(Fy(/)||^ and ||»||^..= X;||(V*)«(»)||^. 

i=0 ~ t=0 

Indeed, V, V* and V, V* share many good properties; we summarize the analog of Lemma 
EES ESI E2I and S3] for F, F* in the following. 

Lemma 3.9. (1) If f G satisfies gm-i G L| /or some rea/ number m, m + A; > |, i/ten 



/ 1,2 <N *7 ||2 

2 s II^tIIl 2 



^4/so, if g satisfies G L| and G L| /or some real number m, m < k + ~, i/ten 

II A||2 < II 1 12 , || g i|2 

1 1 £m I II 2 — II £m-l I Il| 1 1 £m-l I Il| ' 

/// G Z fc ' [fcl+1 and 5 G Y kM+l , then we obtain 

ll^rlUf ^ C||/||— fc.ffci+i and ||^|U°° < C||5lly&.W+i • 

f5j Lei < j < k — \ be a given nonnegative number. If f G x^'^ +1 and g G 
f/ien we obtain 

ll^jlkf < C||/ll^fe,m+i « nd lljjl^f - c 'll5'llyfc.rji+i • 
(^J Product rule for V , V* : 

V*f = -Vf + 2k f - V(fh) = V(f)h + fd£h, V*(gh) = T(g)h - gd^h. 
Now we show that two norms induced by V, V* and V, V* are equivalent. 
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Proposition 3.10. Let 4> be given so that \ \C k( P\ \x k ^ k ~\+ 2 — A and ^ < 4 < C for positive 
constants A, C . Then for any f S X k '^ k ^ and g £ Y k '^ k \ there exists a constant M depending 
only on A, C so that 

^H/ll^.w < \\f\\ 2 xk ,m < M\\f\\^ k , w , ^\\9\\ 2 7 Kw < Ibllkw < M \\9\\y k ,m • ( 2 o) 

Proof. We will only prove the second inequality (<). The other one can be shown in the 
same way. Let us start with V f and V*g. 

1 r 6 2k , 6 2k — r <A 2 *\ 6 2k — t 2k V(£ k cb) 6 2k , 

Vf = ^ d ^f] = §ff*7 + d^]f = ^Vf + {y^^pf- - 2fc^}/ 

Next, note that by the product rule of Lemma 13 . 8 1 and [3 . 9 \ the higher order terms (y) l f and 
(V*) % g can be expanded into the following form in terms of (V) J f and (V ) 3 g for j < i: 

i i—1 
3=0 3=0 

where for s = 1 or 2 

r=0 ^ h+-+l p =j-r 9=1 ^ ^ 

Zi,...,Z p >l 

for some functions C rs , C^.../ pS which may only depend on k, f and 4 and therefore C rs , 
Ci x ...\ s are bounded by some power function of C. In order to show (|2U[) . first we rewrite 
(V)7 and (V*Yg as 

j=0 5 j=0 ^ 



Since 1 1 ^- ^ \\l 2 an d II gj ^ Hl 2 are bounded by ||/||_ fc j and \\g\\y k,j respectively and 



|(Z1!Z/|I T2 and \\Qq£i\ 

moreover, by Lemma I37FI H^^IIl^ is bounded by ||£ fc 0||xfcj+2 for j < \k], by adding all 
the inequalities over i < \k~\ we obtain the desired inequality as well as the desired bound M 
as a function of A and C. □ 



The above equivalence dictates the linear character of higher order energy. Note that we 
cannot deduce the same result for the full energy £ fc >r fc l+ 3 d ue to the nonlinearity. 

4 V, V* a priori energy estimates 

This section is devoted to V, V* energy estimates to get the following a priori estimates, a 
key to construct strong solutions. 
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Proposition 4.1. Suppose 4> and u solve (|1Q|) and (jlip mf/i £ fc, r fe l +3 (0, -u) < oo. If we 

further assume that 

]^<% <C for some C > 1 , (21) 
we obtain the following a priori estimates: 

^-£ k >W+ 3 (cf>, U ) <C(£ k >W+ 3 (<f>,u)) 
dt 

where C(£ k '^ k ~\ +3 ((j),u)) is a continuous function of £ k '^ k ^ +3 (4>,u) andC. Moreover, the a pri- 
ori assumption (121h can be justified: the boundedness of £ k ^ k ~^ +3 ((f), u) imply the boundedness 
<>f j- 

In order to illustrate the idea of of the proof, we start with the simplest case when k = 1 
of which corresponding 7 is 3. In the next subsections, we generalize it to arbitrary k > ^. 

4.1 A priori estimates for k — 1 (7 = 3) 

When k = 1, the Euler equations read as follows: 



(22) 



The operators V and y* take the form: 

1 6 2 6 2 1 

In terms of V and V*, (|22|) can be rewritten as follows: 



3 

The energy functional £ 1,4 (<f>, u) reads as the following 



d t m-V*(Zu)=0 

1 (23) 



4 

^■ 4 (0, «) = y ^i 2 + i^i 2 ^ + flj ^(^)i 2 + icn'(e«)i 2 <« (24) 

Before carrying out the energy estimates, we verify the assumption f) 2 1 1) . In order to do so, 
we examine each term in the energy functional (|24p . Let us start with V, V*V, VV*V, 
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V*VV*V of £0. 



i 



One advantage of V and V* formulation is that L|° control of (ft is cheap to get: since 

J £0 • vmdi = j 4>d^n = \<t>\ 

we deduce that 

'1,4 



~<\mu + wmwu <s^(m (25) 

The boundedness of also follows from the boundedness of £ ' (<^>, u): first note that 

F^ = 3^ 

and then by applying Lemma [3^61 to g = V(^<p) when k = 1, we deduce that ^d^cf) is bounded 
and continuous, and in turn d^<p is in under the assumption pip . In the same way, by 

applying Lemma 13.61 to / = V*V (£(/)) when k = 1, we can derive that d^[^d^(j)] is bounded. 
However, we remark that this does not imply that <9|0 is bounded in our energy space, since 

it is not clear how to control Thus we keep the form as they are rather than try to go 

back to the standard Sobolev space. 

Next we turn to u variable. We list out V*, VV*, V*VV*, VV*VV* of £u. 

6 2 
4 



1 

2 ^ 



VV*^u) = --d^d^u] 



V*VV*{iu) = ^d i [- 2 d i {^u 



(26) 



Note that can be estimated in terms of d^u via the equation ([22 



o f 6 = — "rr 0c u - 



Apply Lemma 13.61 to deduce that vv = -^d^[(j) 2 dt(p] is bounded and continuous if 
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£ 1,4 '(4>, u) is bounded. Letting h be -^d^[4> 2 dt(j)], we can write 4> 2 dt4> = Jq i' 2 hd^, and therefore 
we conclude %^ is also bounded by S 1,4 ^, u) and C. 



Writing |? as 



|(«) = f(o,e)+/'f (r,e* 



we conclude that for a short time, the boundary behavior of | is preserved and in particular, 
this justifies the assumption ([2Tj) . 

We now perform the energy estimates. The zeroth order energy energy is conserved as 
given by Apply V and V* to ([23} and use V t (£4>) = §$V(£^) 

$v(C0)-3W*(£v) = o 

^v*(eu) + iy*y(^) = y t *(?u) (27) 

Multiply by §V(£0) and V*(£u) and integrate to get 

lit J l^ V ^)\ 2 + l^*(^)| 2 ^ = / V?(fr) • V*(£u)dt (28) 
Apply V* and V to ([27]) to get 



2V*V t V*(£u) ■ V*W*($u)d£ + J V t *VV*(£u) ■ V*W*($u)d£ 



(30) 



dtV*V(£<f>) - 3V*VV*(£u) = V t *V(£<j>) 

1 (29) 
d t VV*^u) + ~VV*V(£4>) = VV t *{iu) + V t V*(£u) = 2V t V*(£u) 

Multiply by \V*V(£4>) and VV*(£u) and integrate to get 

\jt J l\ v * v ^)\ 2 + \w*(tu)\ 2 dti = j \v;vm ■ v*v(cm 

+ J 2v t v*{e,u)-vv*{iu)di 

Apply V* and V to ([25} to get 

$v r v*v(£0) - 3yy*yy*(^u) = yy/y(^) + V5V*v(^) = 2V r t v*v(£0) 

1 (31) 
d t V*VV*(£u) + = 2V*V t V*{£u) + F/W*(£«) 

o 

Multiply by \VV*V(£<j)) and and integrate to get 



(32) 
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Apply V* and V to (J3TJ) to get 



d t V*VV*V(£cj)) ~ 3V*VV*VV*(£u) = 2V*V t V*V (£(/)) + V?W*V (£(/>) 
d t VV*VV*(£u) + \vV*VV*V(Z</>) = 2VV*V t V*(£u) + 2V t V*VV*(£u) 

Multiply by \V*VV*V '(£<(>) and yy*l/F*(£u) and integrate to get 

~ f -\V*VV*V(&)\ 2 + \VV*VV*{£u)\ 2 d£ 
2 at J 9 

= J^V*V t V*V(£<f>)-V*VV*V(Z<f>)dS + J ^V t *VV*V (£0) ■ V*VV*V(£<t>)d£ 
+ J 2VV*V t V*(£u)-VV*VV*(£u)d£ + J 2V t V*VV*(£u) ■ VV*VV*(£u)d£ 



(33) 



(34) 



In order to prove Proposition 14.11 it now remains to estimate the following nonlinear terms 
coming from the energy estimates in terms of the energy functional 



2V t V*(£u) • VV*(£u)d£, J ^V t V*Vm • VV*V(£(j>)d£, 

J 2V*V t V*(£u) ■ V*VV*(£u)d£, J -V*VtV*V{£4>) ■ V*VV*V{£(f>)d£, (35) 

J 2VV*V t V*(£u)-VV*VV*(£u)d£, J 2V t V*VV*(£u) ■ VV*VV*(£u)d£ 

Our goal is to control these terms by our energy functional (|24p . All of them contain dt<j> 
and its derivatives with suitable weights. The estimates of ft</> related terms can be obtained 
through the equation ([22]) by estimating V*, VV*, V*VV*, VV*VV* of £u in terms of d t <j>. 
Before going any further, let us try to get the better understanding of them. First, we rewrite 
V*, VV*, V*VV*, VV*VV* of£u, namely ®, in terms of 



V*{£u) 



£ 



■d^u = £d t 



1 „ r ,0„ „ (tP „ rft<, 



VV*(£u) = -diWM] = +3^^ 

£ £ £ 



£ 4> 



+ vm d 4 



V*VV*(£u) 



j d d-^ d d— ]] " 3 y<%[^ " 3 

1 „ r (j) 6 „ ,$0,. o 2 _ r (/> 2 _ ,,ft0 



/> 2 </> 2 



ft. 



V 9 * 

^2" 



ft 



ft< 



ft 



(36) 



Of 



£ £ 
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As we can see in the above, the term ^ and its derivatives naturally appear and there are 
many ways to write them. The key idea is not to separate them randomly when distributing 
spatial derivatives, but to find the right form of each term. The boxed terms in (|36p have 
been chosen in such a way that the remaining terms in the right hand sides have the better 
or the same integrability as the left hand sides. 

We analyze the most intriguing full derivative terms V*V t V*V (£</>) and VV*V t V*{^u). 
Other terms can be handled in a rather direct way. 

Claim 4.2. 

\\v*v t v*vm\\q <Ci(£ 1,4 (<M) 

where Ci(£ 1,4 ((f), u)) is a continuous function of S 1 ' (<j),u). 
Proof. 



1 d t c/> ^ J> 2 

EE (/) + (//) +(///) 

Since (J) = ^F*W*F(£<£), (J) is controllable: 



L|<ll-^lkrll FWV (^)Hii 



The second term is written as 



From (|36|) . we get 



(//) " "3T^ [ — 1 c — 



J dd ~ ] - —* G L? • 



And since 1 1 VV 1 1^ 2 < ||V*VV*V(^)|| £( 2, (I/) is also controlled by the energy. Now we 
turn into (III)- First, we note that by Lemma 13.31 and Lemma 13.61 

-^4 d ^\ e L\ and d^d^} e Lf. 
Next let us take a look at the other factor and rewrite it by using boxed terms in (|36p: 

4P 1 5*0 6 1 6 ^ 1 J 6 dt<f> J? 

■^ d d^[—]] = —%[——%[_]] = ^dd^[—}]-6^d^-^[—] 

= W (*0 + — 1 _ + 2-| — | — 
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Thus (III) can be rewritten as follows: 



</> V*VV*(tu) <j> 2 V{0) VV*(Zu) 1 <j? 
J 1 d d£l d ^} "2— 1 |^£[|2 

— + a^l—l — 



Hence can be controlled by £ l,i (<p,u). □ 

Now let us move onto VV*VtV* (£u) . The treatment of this term contains another flavor. 
Claim 4.3. 

\\VV*V t V\iu)\\l l <e 2 {8 l '\cl ) ,u)) 
where C2(£ 1,4l (4>,u)) is a continuous function of £ 1 ' i ((j),u). 

Proof. We use the continuity equation in ([22]) first to deal with VtV*(£u). Since V*(£u) = 
£d t <l) and V{£d t <l>) = |d ? [c/> 2 <9 t </>] = \[(j) 2 d^d t (j) + 2^d^d t c/)], we can write V t V*(£u) as the 
following: 

\v t v*{iu) = -M\d t <p\ 2 ] = -Ad^dt^ 2 + 2<pd t d^d t <p] = 2^vm<t>) - 3^|^ 

2 ? ? 9 ? 

Apply F*: 

2 c <p ? r 

V v ' 

(*) 

Note that (*) reduces to: 



1 «u3^i odttP^ 2 *! ^aA&i 



Apply V: 



±W*V t V*(tu) 

= ^rvrN - 2^|^]|-] + s^^llf I 2 ] 

= (I) + (//) + (III) 
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We rewrite (I), (II), (III) as follows: 



(in) = 3 d^d^d^}]\— | + e- ^ • — 



It is easy to see that (/) and (III) can be controlled by the energy functional. On the other 



hand, in order to take care of (II), a special attention is needed. Since — V K ^ is bounded 



by VV*VV*(iu) in L|, we obtain 



1 „ r ctP„ r d W v ^ T 2 



h =-fAl^d^)]eLj. (37) 

Thus the second term in (II) is bounded by the energy functional. To prove that the first 
term is also bounded, we claim 



By using ([37]) . rewrite |t<%[^] as follows: 

] = / c hdC 

Applying Holder inequality, we observe that 



o 



\^[^]\ 2 <e\\h\\i 

Hence we get 



This finishes the a priori estimates for k = 1 as well as the claim. □ 



4.2 The case when ~ < k < 1 (7 > 3) 

In this subsection, we prove Proposition 14.11 for the case | < k < 1 and s = 4. First, by 
Lemma 13.61 one finds that d^(f>, 4, and ^ are bounded by the energy functional £ k,i (<p,u). 
We recall the reformulated Euler equations (|1U|) , One can apply V, V* alternatingly as in the 
case k = 1, and integrate to get 

4 

/ ( 2 fc+ ip lOO^)! 2 + im*(^)| 2 ^ < mixed nonlinear terms, 

i=i 

where the mixed terms have the same form in (|35p . Thus it suffices to show that those 
nonlinear terms are bounded by £ k,4 ((f),u). We focus on the intriguing term VV*VtV* (£, k u) 
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as well as V*VtV*(£ k u), VtV*{^ k u). As we saw in the case k = 1, in order to treat the 
mixed terms, the careful analysis of dt4> terms is required. First, take a look at V*(£ k u) and 
VV*{i k u). 



V*(£ k u) 



x2k 

1 - F^l J. (01? 4. ^\ M 



2fc+1 /H^l 



Because 



and 1 1 — — -p— ^ - 1 1 L 2 are bounded by £ fc ' 4 as an application of Lemma 



13.31 we deduce that ||</> fc c?g[^]|| L 2 is also bounded by £ fc ' 4 . And by Lemma [3T6l we obtain 
(j)d(\p^-\ € L|°. Now let us compute V(V*(£ fc u) and write in terms of the energy: 



^^«*«) = ^ 2fc+1 l7l 2 ] = "-fi^^-f] + (2* + 1%^} 



1 



It is clear that this term is bounded by the energy functional. Next we write out V*VV*{t} 



= -^^^-^plp^f 



(38) 



1 „> 4fc+2 £> t <£ 



Thus we note that the boxed term ^jd^[ - d^[^-]] is the right form of the second derivative 
of dt<j> of which structure we do not want to destroy. Here is V*VtV*(^ k u): 



}] 



^ i „ r <A 4fc+2 



(» + i)l^-^[pM 



It is easy to see that the first and third terms are bounded by the energy functional. For the 
second term, writing it as 



</> 4fc+1 dtd> 2 _ ^ k+l At k 
2 -^-\ d d—\\ -- 2 ^3jr d d—\^ d d—\ 
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we deduce that it is also controlled by £ fc ' 4 , and in result, we conclude that V*VtV*(£ k u) is 
bounded by £ M . Next VV*VV*(£ k u): 



i (p- 1 <* 4fc+2 



Note that the boxed term is bounded in L| by the energy functional. Now we write 
VV*VtV*(^ k u) in terms of the boxed terms as well as the energy: 

—VV V t V (C n) = -2— • —^[-^-^[-^fl^—]]] 



-(*) 



-(2k + 1)|*4|> • ^[Ofel&Wl - (4b + 2)^^^ ■ ^ 



It is clear that except (*), the first factor of each term in the right hand side is bounded in 
and the second factor is bounded in L| by £ fc ' 4 . In order to show that (*) is bounded in L|, 
first note that by Lemma I3T31 1 1 v vv ^ ") 1 1 £2 anc l 1 1 v v j€ ^ \\ L 2 are bounded by £ fc,4 (0, u) 
due to the fact k > i, from the relation (|38|) . thus we get 



and in turn we obtain 



4fc+2 



||^|^[M]|2||2 2 < | N |4 2 /' 1 ^ + (4 fc+ 10)-(8 fc+ 8)^ 

« Uo 



Note that the last integral is bounded. Therefore we conclude that | |yy*Vty*(£ fc u)| \ L 2 is 

bounded by £ k,i . Similarly, we deduce the same conclusion for other mixed terms and it 
finishes the proof of Proposition 14 . 1 1 for | < k < 1. 

4.3 The case when \k] > 2 (1 < 7 < 3) 

We now turn into the general k. The spirit is the same as the case when k = 1: we need to 
carry out estimates and nonlinear estimates. For large k, however, the number of mixed 
nonlinear terms increases accordingly and it is not an easy task to work term by term. We 
will present a systematic way of treating those terms involving derivatives of dt<f>. 

The following lemma is a direct result from Lemma 13.61 and [ 3.71 an d it is useful to justify 
the assumption (|2~T1) in S k '^ k ^ +3 ((j),u). 
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Lemma 4.4. (1) We obtain the following L|° estimates: 

where C-z(£ k ^ k ^ +3 (4>,u)) is a continuous function o/||^||x,«> and £ k ^ k ~^ +i ((j),u). 
(2) For 3 < i < \k] + 1 

ii W(^) |i N (n!(iH| < r „^ 

where C^(£ k ^ k ^ +3 (4),u)) is a continuous function o/||^||i|° and £ k, ^ k ^ +3 ((j),u). 

We start with L|° estimate of <j). We list out V, V*V, VV*V of for references. 

• V(£ k <l>) = ^d^ 2k+1 ] = (2k + l^dtf 

cf? k ^ 



^) = -(2Hl)^t% 



• VV*V(^) = -(2k + l) Tk d,[^-d^d^)] 
One advantage of V and V* formulation is that L|° control of eft is cheap to get: 

Applying Lemma 14,41 we obtain that d^cj) is bounded and continuous if the assumption (|2ip 
holds and £ fc, ^ +3 (</>, u) is bounded. 
Next we turn to u variable. 

• v%eu) = -^u 

1 rh 4k 

. Wfc*ti) = --d^dsu] (39) 
Note that ^ can be estimated in terms of d^u via the equation ([5j): 

j,2fc \l*(ck 



£2fc « £fc 

By Lemma 14.41 we deduce that — = -^d^[(f) 2k dt(f)] is bounded and continuous if 

£ k '^ k ^ +3 ((f>,u) is bounded. Letting h be ^<9£[0 2A: <%(?!>], we can write 4> 2k dt4> = Jq £,' 2k hd£,, 
and therefore we conclude is also bounded. By the same continuity argument as in k = 1 



case, we can verify the same boundary behavior for a short time and the assumption (|2ip . 
We now perform the energy estimates. From (|l(jp . we have the conservation of the zeroth 
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order energy: 

i d r 1 



2 y 2k + 1 1 

Apply V and to {TO} and use V t (£ k <f>) = jE&dtV {£ k <}>) to get 



d t V{£*<t>) - (2k + l)VV*(i k u) = 



9 * y * ( ^ n) + 2TTi FV( ^ ) = Vt *^ u) (40) 

Apply V* and V to (gDJ) and use FF t * = I^V* to get 

By keeping taking V and V* alternatingly, one obtains higher order equations for (V) l (^ k (j)) 
and (V*Y(^ k u) for any i. Indeed, the mixed terms in the right hand sides can be written in 
a systematic way. For i = 2j + 1 where j ' > 1: 

j-l 

^(F) 2i+1 (^V) - (2ife + l){V*) 2j+2 {i k u) =2 s ^{V*) 2j ~ 2l ~ 2 VtV*{V) 2l+1 {t k <t ) ) 

1=0 
3-1 



(42) 



^(y*) 2j+1 (e fc ^) + ^^(^) 2j+2 (cV) = 2y(yfj- 2| -v 1 f(f) 2! K 

2k + 1 ^-^ 

1=0 

+ v t *(y*) 2j (t k u) 

For i = 2j where j > 2: 

3t(V0 2j '(£V) - (2Jfe + l)(y*) 2j+1 (e fc n) = 2Y,(V*) 2j ~ 2l ~ 3 V t V*(V) 2l+1 (Z k (i ) ) 

1=0 

+ V t *{V)V-\e$) (43) 
1 i=l 

d t (v*) 2j (t k u) + — — (^) 2i+1 (eV) = 2V(r)^^V ( r(y*) 2, (^) 

zfc + 1 ^-^ 

1=0 

Our main goal is to estimate the mixed terms in the right hand sides of (|42p and (|43p . Note 
that the most intriguing cases seem to be when 1 = 0, where the most spatial derivatives 
of ^ are present. Before getting into the estimates, let us get the better understanding of 
the effect of the operator Vf. First, recall (V*)*(£ fc u)'s. They give rise to dfcj) terms via the 
equations ([5]) and furthermore they predict the right form of spatial derivatives of dt<j) terms. 
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We borrow the computations from the previous section for i < 4. 



fh 2k 

. V*(en) = -^d^u = ed t( t> 



» 2fe+1 Acf> ,<f> 2k A<f> 



IT (44) 



J 2 ** 1 ^ 



Now let us compute VtV*(£ fc u) and write in terms of the energy: 



=2 d 4w*{e U )-\ d 4\ 2 v{z 

6 4> 



Next we write out 

^0 



Thus we note that |^9g[ ^ 2 fc" ^[^]] i s the right form of the second derivative of that 
we do not want to destroy. Here is V*V f V*(£ k u). 



do not want to destroy. Here is V*V t V*(£ k u). 

= -f>^<]] - **p<ir - <* + u£*[£*<f- <46> 



-}] 



Next ww* 



1 J,2ifc-1 j,4fe+2 o~T q j, <2fe Q i 



vv*v(e<p)?f+pv*v(e<t>)-d^] 



(47) 
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In turn, VV*V t V*(£ k u): 

i k M 1 2fc -i 4fe + 2 a</> 

-vv v t v a u) = -2— ]]] 

" + Dlf I 2 • ^8^3^ - (4, + 2)^9^ • ^] 

As in the case n < k < 1, we chose the above specific expansion of V*VV*{^ k u) in order 
to find the right expression of the second derivative of dt<f> and then estimate the mixed 
term V*VtV*(£ k u). This is because the nonlinearity of ((5J is prevalent at this level in V, V* 
formulation and an arbitrary expansion may destroy the structure of ©. For higher order 
terms, we employ a rather crude expansion in a systematic way. In the below, we present 
a representation for (V*Y(^ k u), get the information of spatial derivatives of dt<p, which we 
denote by Tj's, and derive the estimates of (V*y~ 2 VtV*(^ k u). 



Representation of (V*y((, k u) and T; for i > 2: First, we define Tj for 2 < i <= \k] + 3: 

n = '^8 tl f ]; T 3 = -^[^[f ]]; 7S = W «T 3 for « > 4 (48, 
We chose T 2 and T 3 so that 

(V*) 2 (t k u) = T 2 + (OV*) = T 3 + (y) 2 (eV)^ ■ 



In other words, they are boxed terms in the above as well as in the previous sections. Note 
that by Lemma 14.41 

^€Lf. (49) 

The validity of Tj will follow from the construction of T, in the below. Now we claim for any 
i > 4, (V*)*(£*«) has the following representation which extends the case i = 2, 3: 

2 — 2 

(y*r(e fc n) =Ti + (vy- 1 ^^ + I>-/z) ( 5 °) 



where 



i-j-2 



r=0 ^ iiH Yl v =i~j~r q=l * V 

li,...,l p >l 

for some functions C r , C^...^ which may only depend on | and | and therefore C r , Ci v ..i p 
are bounded by ||^||l°° and Furthermore, T{ and the last term in ([50]) have the 

following property: for each 4 < i < \k~\ + 3, 

Ib^E*^!^ ^d are bounded by S^+^u) , (52) 
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and for 4 < i < \k] + 1, 

i— 2 

ibsjW£*-i r ili*r and Wjrn^^r are bounded h * ^ ,W+3 (<M- w 

? 3=2 ? 

In particular, each Tj is well-defined. We are now ready to prove the representation formula 
(j50|) . From the previous computation ([17)1 . 

(n 4 (e fc u) = r 4 + ew*)^ + (y) ^ V) r 2 . (54) 

Setting $4_2 = ^gil ^ , the formula ()5Q[) holds for i = 4. Moreover, by Lemma 14.41 the 
properties (f52j) and (f53l) are satisfied for i = 4. For i > 5, based on the induction on z, we 
will show that V or V* of each term in (|50p can be decomposed in the same fashion. With 
the aid of Lemma l3.8l applying V or V* of the second term iy) l ~ l {^ k (f))^- in the right hand 
side of ([50]) . we obtain 



Note that the right hand sides are of the right form. In the same spirit, we can apply V or V* 
to the last term in (i50|) . Here is the estimate of the simplest case when j = 2, r = 0, l\ = i—2. 
Consider i = 2j + 2 or i = 2j + 3. 

• M ?fc7 i 2J = - n. , -i Tkl Tkl ^2 TT-7 J2 + 



£ k 4> ' 2fc + l £ fc </> £ fc </> £ fc </> 

Each term in the right hand sides has a desirable form. From Lemma 13.81 we deduce that 
when V or V* act on a function h, depending on that function, they can yield | or V ^ k ^ ^ 
or V/i or Therefore, V or V* of other cases of &i-jTj falls into the right form of the 

case when i + 1. Next we verify (|52p and (|53p . Let s > 4 be given. Assume that they hold 
for i < s and we first claim that 



j (*+l)-2 

L W+3 -(s + i) E ^ bounded by £ k >W +3 (</>,u) 



This can be justified by counting derivatives and | and distributing | factors in a right way. 
The highest derivative term with appropriate factor of | will be taken as the main L| term 
and others will be bounded by taking the sup. Let z = max{j, l\, . . . , l p } for each term. If 
z = j, we consider ^l 3 _ j whose L| norm is bounded by £ k '^ k ^ +3 ((f), u) from the induction 
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hypothesis, since j < s — 1, as an L 2 term. Now it remains to show that the rest of factors 



£[k]+3-(s+l) L-t pr 11 tkx 

S r=0 s Ii H hip= — i — r 9=1 

Zi,...,Z P >l 

are bounded. It is enough to look at the following: for each r < s — j — 1 

ii+--+/ p =s+l-i-r?=l S Y 
h,...,lp>l 

which is clearly bounded by £ k '^ k ~^ +3 ((j), u) because of Lemma 14.41 and since k < \k~\. When 
z is one of L's, one can derive the same conclusion. Thus from (|50p . we deduce that 



[ | g-ffc l +3— (s+i) \ \l 2 ^ s a ^ so bounded and it finishes the verification of ([52]) . For the bound- 
edness for s + 1 < \k~\ + 1 

1 (-+1J-2 

IL rj fc1 +2 _ (s+ i) S *(»+i)-i T illir 
^ i=2 

we employ the same counting argument except that we have to use the fact , ^ fc1 + 2 -j £ 
for 3 < j < s. This finishes the induction argument. 

These Tj's and their properties will be useful to estimate the nonlinear terms involving 
VtV* . Based on the representation formula (f50l) and Tj's, we further study the representation 
of more general mixed terms (V*) l Vtf for i < [~£f| +1. First Vtf can be written as 

Note that the right hand side has the same structure as the last two terms of the right hand 
side of ([Ml), in letting / be (V) 2 <j>) . Thus we can apply the same technique to obtain the 
following: for each i < \k] + 1 

-(vyv t f = J^(vy+ 1 f + ^( m+ 2)^mvr m f. (55) 

^ m=0 j"=2 

Proposition 4.5. (Nonlinear estimates) The right hand sides in (|42j) and (|43|) are bounded 
in L 2 by a continuous function of £ k '^ k ^ +3 (<f), u) and C. 

Proof. We only treat the most intriguing term (V*) l VtV*(£ fe u) for i < \k~\ + 1. The idea is 
to estimate it in terms of Tj's and to use their properties. First, note that V*VtV*{^ k u) can 
be written in terms of T 2 , T3 as in 



— V*V t V*(£ k u) = 2^T 3 -2-t- ,rr ' 12 1 ^*Tr^fc^i^,> 
2k y ' (p £ k - 



—V*V t V*(Cu) = 2^-T 3 - 2 7r7 |T 2 | i + V*V(C(f))\-^-\ 2 (56) 
We would like to compute {V) % V*V t V*(^ k u) for 1 < i < \k]. Consider i = 1. We apply V to 
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each term in (|56p . Since 

d t 4> £ k 

we have 

V(— |T 2 I 2 ) = ^T 2 • T 3 - **±* J_^M W 

y(y*y(e^)|^| 2 ) = yy*y(£V)l^l 2 + • v*v(e^ 

Thus we deduce that VV*VtV*(£ k u) is bounded in L| by the energy functional. In the same 
vein, by keeping applying V* and V, for any i > 2, we can write (V) l V*VtV*(£ fe it) as follows: 

irymr^u) = 2^t, +3 + (vr +2 (^)l^l 2 + ^ -E*^-^ 

j=2 (57) 

i+2 i+2 s-2 . v ' 



where $ is given as in (|5ip with possibly different coefficient functions and Cj and C s j are 
some functions bounded by 1 1 1 1 1 l|° and 1 1 J | . This formula ([57]) can be also obtained by 
plugging ([3D]) into ([35]) . 

Now we claim that for each 1 < i < \k],\\(VyV*VtV*^ k u)\\ L 2 is bounded by £ fc 'W+ 3 ((/>,u). 

The first three terms in the right hand side of (j57|) are bounded since they are of the form as 
in (|50p multiplied with The rest terms consist of quadratic or higher of energy terms or 
Tj's. For each term, the highest derivative with appropriate factor of | is considered the main 
L| term and other factors are bounded by taking the sup. This can be done by employing 
the counting and distributing | argument as well as the estimates of Tj'a ([52]) and ([53]) as 
before. □ 



5 Approximate Scheme 

In this section, we implement the linear approximate scheme and prove that the linear system 
is well-posed in some energy space. 

Let the initial data <^>o(£) and ^o(0 of the Euler equations (J3J) be given such that < 
^ < C for a constant C > 1, and £ k ^ k ^ +3 (4> , u ) < yl for a constant A > 0. Here 
5 fe <r fe l+3(^, 0) Uq ) i s the energy functional ([21]) induced by </>o- Note that from the energy 
bound, we obtain G We will construct approximate solutions <fi n (t,0 and u n (t,£) 
for each n by induction satisfying the following properties: 

4>n\t=o = <fio, u n \ t =o = u ; 4>n\z=o = «n|^=i = 0; < < C n , for C„ > 1 (58) 
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Note that (fto, u o automatically satisfy ([58]) . Define the operators V n and V* as follows: 
In addition, we define commutator operators {V n )t and (V*)t- 



VM=eiW%rf)> K(9) = -^ds(- k g) (59) 



We define 9t^>o through the equation by 

dt<f>o = -|§f^ u o- 

For the linear iteration scheme, we approximate (V) 3 (^ k 4>) = G and (V*) 3 (^ k u) = F and the 
equations (|42p when j = 1, instead of <j> and u themselves. Let 



For each n > 0, consider the following approximate equations 



8 t G n+l - (2k + l)V n F n+1 = 4 

v c 

2k + 1 n 



dtF n+ i + — — —V*G n+ i — J, 



where and are given as follows: 

4>n Z k 4>n 4 

The initial and boundary conditions are inherited from the original system: 

G n+1 | t=0 = WVo(C'Vo), F n+1 \ t=0 = V *VoV *^ k u ), G n+ ik=i = 0. 

Note that F n+ \\^ = o is built in the equations (|60|) , In turn, from these -F n+ i, G n +i, we define 
Ai+1> H n+ i, 4> n +i, u n+ i: 



JO ?n+l Jl 



Note that we have used the boundary condition at £ = 1 in order to invert V*. Also note 
that from the above definitions the following identities hold 

V n V*D n+ i = G n+ i, V*V n H n+ i = F n+ i, D n+ \ = V n+ \(^ k <p n +i) , H n+ \ = V* +1 (^ k u n+ i) . 

In view of Proposition 15-H it is easy to deduce that (V*) l D n+ i and (y n ) % H n+ x for < i < 
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\k~\ + 2 are well-defined, namely bounded in L|. Also V n +i and V* +1 are denned as in ([59]) 

with (f) n +i- The right hand sides J*, of ([60j) are approximations of 2(V n )tV*V n (£ k 4> n ) and 
2F n *(y n ) t y n *(e fc n n ) + {V*) t V n V*^ k u n ) in the following manner: 

^(K)^K(CVn) = -(2* + l)^d d ^^^d^ W ] 

1 . .ft^ <j% _ ^ r ^ ni TCK(e fe 0n) 



= - (2* + 1)|5^["^ • ^[pF^]] = -^fKTCK(e^n) + ^[-^ ,, 
, 4> 2 n a { d t ct) n V*D n 

^(WJ^n)~-2^^[^^]]-2^|^[ 1 -]| +|_| ^A, 

2/C n 0„ 

In particular, note that the approximation of V^(Ki)iV^(£ fc u n ), which has the strongest 
nonlinearity in view of the a priori estimates, is based on the expression ([46 p . Also note that 
the equation (|60jl converges to ([42 p for j = 1 in the formal limit. 

We define the approximate energy functional at the n-th step: 

^ r l 

i=0 J i^ + lj ( 61 ) 



1 



"Gn+il P fcjfc] + | l-^n+1 1 |y fc.ffel 



2/c + 1 ~ 1 y„ ' "-'^"jc- 

where A^' s and y n M denote and y fe ' s induced by V„ and V* in ©. 

We now state and prove that the approximate system (j60p are well-posed in the energy 
space generated by X n , Y n under the following induction hypotheses: 

(HP1) £ k < oo and <j) n , u n satisfy 



(HP2) when k < 1, ||%^|| LOO and ||£3J^]|| L oo are bounded by B, 2r_ i; 

V, fc^l ll^ll ll^ll || T i,n „ nWCeVn),, 
whenA:> lJ|^|| ir J|_|| Lr J|___|| ir) || ___ || L?) 

for 3 < i < [fc] + 1 are bounded by S^-i, 
(HP3) ^ an d ^ in TO are bounded in Y" n fc ' W , A^ fc l respectively. 

Here Tj jTl is the Tj defined in ([48h where <^> is taken as <^ n . 

Proposition 5.1. (Well-posedness of approximate system and regularity) Under the hypothe- 
ses (HP1), (HP2), and (HPS), the linear system (j60p admits a unique solution (G n +i,H n +i) 
in Yn' ^ , Xn ^ space. Furthermore, we obtain the following energy bounds. 

Jo 
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where Cs(£'^_ 1 , £^+i) * s a continuous function of S^-Ij ^n+i an< ^ Go- 

Proposition 15.11 directly follows from Proposition 7.1. In the next subsection, we verify 
the induction hypotheses. 

5.1 Induction procedure 

In order to finish the induction procedure of approximate schemes, it now remains to verify 
the induction hypotheses (HP1), (HP2), and (HP3) for n + 1 as described in Proposition 15.11 
The spirit is the same as in the a priori estimates. However, since n and n + 1 are mingled in 
the energy functional (I6ip . the verification of the induction hypotheses needs an attention. 
We only treat the case when k = 1. Other cases can be estimated similarly. First, it is easy 
to see that </> n +i and u n+ \ constructed in the above satisfy the initial boundary conditions in 
(|58p for n+1. The boundedness of &n±l will follow from the continuity argument by using 

the estimate of 9t ^ +1 , Write the equation for <9 4 ((^ +1 ) from the definition of (f> n+ is: 
3^ +1 ^ n+1 = 4 /V) 2 T P tiG^idtidbtf 



Jl Yn JO 

(O 2 / ? 4 [^\idtG n+1 d^ 2 de 



r I x 2 / ^2_ 

n Jl 4>n JO 

^ - re $ 0t * /* 



JO Jl <Pn Jo 



(62) 



3 / ^') 2 / 5^+1^2^' - / (O 2 / % / ClJ^lrf 

JO <Pn JO Jl <Pn JO 



The first term can be controlled as follows: since 

. A„G n+l 



^iG n+1 d^<^\\^\\ L 2 (63) 



it follows that 



(O 2 / ^SP/ eiG^.ideici&de' < € 

Jl 0n JO 



2,| £ ,|4 ndtffinn n G n+ i 

2 rm 7",oo r 2 



Thus ^ Jq(C') 2 Ji ^ 2 ^ n Jq 2 iiGn+id-iid^di' is bounded by £ n +i and the previous energies. 
For the second term, note that 



f & Z? ^ II £ l|2 I, 

Jl Vn 9n « 



Hence, we obtain 



i / (O 2 / ||^ + i«'<iif i 

4 d Jo Jl 0n 0n 



r oc 1 1 1 1 j 2 . 

L f II £ IIL e 
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Finally, since 

/ 64^6 < £ II— — l|° K- i? 

7o *n e f « (64) 

,,||>n||2 11 0n 1,2 iitfl r a *^i|i n ^n-l^ ii 

the last term is also bounded and therefore we conclude that dt( g^ 1 ) is bounded by S n+ \ 
and the previous energies. Note that the nontrivial contribution near the boundary £ = is 
coming from the second term F n+1 . Since 

-^3 (*) = -|s (°) + y o ) dr = |f (°) + y o ) dr . 



we get 



^ -T||ft(^)|k r < ^ < C$ + T\\d t (^l)\\ Lr 



and in result, we also obtain C n+ \ in ([58]) depending on T and energy bounds. 

We move onto (HP2) . Since = ^f^- ^ , we immediately deduce that 1 1 1 1 mo 

is bounded. Take <9t of 



0n+l 7l <Pn 



^ rftf * (65) 



Therefore the boundedness of | L ~ follows by the same reasoning. 

For (HP3), we now show that J^ + i and J^+i of (|60p at the next step are bounded in 
Y n+n X n+i accordingly. 

Claim 5.2. V* +1 J^ +1 and V n+ iJ% +1 are bounded in L|. 



Proof. The sprit of the proof is same as in the nonlinear estimates of the a priori estimates. 
We present the detailed computation for J^ +1 , which is more complicated than J* +1 . We 
start with V n+1 (^^F n+1 ). 



Vn+ll— 7 ^n+lj = 7^ — h n+ i\ = -zO^—^ T*n+l\ 

^n+1 d t (f> n +l -w r . <An+l t o r^f/Wli 



n+1 y v ' ? - £ 

-v ' L j « „ '^-v- 

dt6 n 4--i F n 



+ 2(2^^1 " 

* v '^v^ 

roc r 2 
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Thus \\V n+i (^±F n+1 )\\ L 2 is bounded by £ n+1 . Next we compute V n+1 ( | ^tl 1 2 F n * +1 D n+1 ) . 



gf0 n+ l 2 1 r ^+l fl r -gn+l ii ^n+l ,21 
VVi+H — 7 1 - ~T d d— 72~ d d 7 J I ~T" 1 J 

4 o i j,4 



Pn Pn+l ' 5T n 0n+l <?Wl 



+ 4( w ^x- w ^)| — I — 



Lf L| 



Thus K + i(|^±i| 2 y n * +1 D n+1 ) is bounded in L|. In order to take care of the rest of J 2 +1 , 
first we claim that 



1 77 [ ^ n+1 ]] 1 1 L 2 is bounded by £ n+1 , £ n . (66) 

r r 4>n+i 



In order to do so, multiply (|65|) by and take to get 
d^ds[—]] + 3d^ n+1 ]^W n+1 — 

A2 a A A2 



d t 4>n+l . „, 3 P'n+lr, , o Ap'n+li 



-(*) 



= ct> 2 n+1 ds<f> n+1 {4 J Adt'-sj Bdi'- J CdH'} + ^{4A-3B-C} 
where A, B, C are defined so that ([65]) can be written as 

+ Wn+ldrfn+ldttn+l = £ 2 {4 J Ad? - 3 J Bd? - J Cd?} . 

The term (V) becomes 

d> 2 1 /■« /■« /•« 

(*) = -^^^+i{-90„ +1 5^„ + i9^ n+ i+c 2 (4 y a^'-3 y Bd t'-j Cd ^ } 

Hence by canceling terms, we obtain 

^[^2~ 9 a^ ]] = - 3 ^n+l^[-T2-O$0n+l]-T + —5— {4A - 35 - C} 



2 * 3 ( 6? ) 

n+172" 1 + -^-i AA " 35 - 6} 



Indeed, A,B,C were treated during the estimates of || ^p±l ||xg°- From the same analysis, 
we deduce that L| norms of ^4, 5, C are bounded by £ n +i, £ n . Thus the claim (|6"6"j) follows. 
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Now we consider F n+1 (^i|%[^f]| 2 ) . 

By ([66]) . the first term in the right hand side is controllable, and again due to ([66|) . since 
l%^[^f]| 2 < e\\pd^dd 9 -^r]]\\l r we obtain 



*.6 



H i ikApn < ll^^^]]!^ 

t, 1 i <?n+l « § C 0n+l « 

and this completes the estimate. For the last term in J 2 +1 , note that 

n+1 n+ i £0„+l ? £ 2 5 n +l 

Thus it remains to show that §<%[^%[^<%[^Jf-]]] is bounded by £ n+ i, £ n . From (1571) . 
one can write it as 



\d^d^d,[^]]] = ±fl^ • fv:D n+1 ] + ^[^(AA -SB- C)] 

= (I) + (II) 

(I) can be decomposed as follows: 



hence it is bounded. For (//), we put A,B,C back into the expression: 

£ <Pn <Pn Jo ? JO 

?n ?n Pn <T Jo 

Therefore, by the estimates (f63|) and ([Mj) of J* ^iG n+ i(i^i and f qiJnd£ii we conclude that 
(77) is also bounded by £ n +i, £ n . This finishes the proof of Claim. □ 
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6 Proof of Theorem 12.11 



In order to prove Theorem 12.11 it now remains to show that <p n , u n converge, the limit 
functions solve Euler equations ([5]), and they are unique. 

6.1 Existence 

First, by applying Gronwall inequality to the energy inequality in Proposition 15. 1\ we can 
deduce the following claim. 

Claim 6.1. Suppose that the initial data (/>o(£) and ito(£) of the Euler equations © be given 
such that ^ < < Co for a constant Co > 1, and £ k, ^ k ~\ +3 ((f)Q,uo) < A for a constant 

A > 0. Then there exists T > such that if for m < n, £^ < |A for t <T, then £ k l+ i < \A 
for t < T and in addition, for all n, < ^ < 2Co for t <T. 

Thus we get the uniform bound of £ k, s as well as the uniform bound of the upper and 
lower bounds of Since the approximate energy functionals (|6ip depend on not only the 
approximate functions G n+ i, F n+ \ but also <p n , the corresponding Banach space changes at 
every step. In order to take the limit, it is desirable to have the fixed space where G n +i, F n+ \ 
live. The plan is as follows: by making use of Proposition 13.101 we prove the equivalence 
between our energy functionals and the energy functional induced by (|18j) so that approximate 
functions have the uniform energy bounds in the Banach space induced by (|18p and thus we 
can apply the fixed point theorem. We present the detailed analysis for k = 1. 

Recall the homogeneous operators V and V : 

and we define the corresponding energy functional £ n +i- 

1 n - - 

£ n+1 (t) = Y,J g\(VjG n+1 \ 2 + \(VYF n+1 \ 2 dC 

i=0 

We claim that £ n +i and £ n +i are equivalent in some sense. We compute £ n +i- V n and V* 
are written in terms of V and V . 

• V*G n +i = -pyV G n+ i 



V n F n+ i = ^VF n+1 + 2(^d^ n - ^Y- 



n+1 



L2 a a2 p 

Note that from the definition of D n 



-jc/^^ni - | — — — | < \\^-d^4> n -i\\ L f + — ^ 2 C n _ 1 ||G rl || i 2 < I + C n „ 1 ||G n || L 2 

where / = |||t^0o||l^ and C n _i is the bound of zn^l as i n (j58j) . Therefore, we deduce that 

< (1 + Af„)F n+1 
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where M n = 5C 2 {C 2 + I 2 + C^„ 1( S n }. To show the the converse, namely £ n +i is bounded 
by £ n +i, we rewrite V and V* in terms of V n and V^. 

• V G n+ i 

• = 

Thus we reach the same conclusion: 

£n+l < (1 + ^nj^n+l 

Note that M n 's have the uniform bound over t < T by Claim 16.11 Therefore, there exists 
a sequence so that G ni , F n; , <fi ni , u ni converge strongly to some G, F, 0, u. Due to the 
uniform energy bound, we also conclude that G, F solve the equations (|42j) for j = 1 and 
moreover (f>, u solve Euler equations ([5]) with the desired properties. 

Next, we turn to the general case. Back to the approximate system (|60p . we define the 
corresponding homogeneous energy functional: 

^n+iW = II 2 £ T ^n+l|l^fc,r*i + ll-Pn+lll^.W ( 69 ) 

From Proposition 13.101 one can derive the equivalence of the associated energy functional 
(|69p and the original approximate functional (|6ip : There exists M n > only depending on 
initial data, C n , C n -\, and £ n such that 

Y^Wjn+l ^ £n + l < (1 + M„)^ +1 . 

By the same reasoning as the case when k = 1, the existence of G, F, 0, u follows. 
6.2 Uniqueness 

In order to prove Theorem 12. H it only remains to prove the uniqueness. Let (<pi,ui) and 
(4>2, u 2) be two regular solutions to ([5]) with the same initial boundary conditions with 

£ k ^ +z (<f> x ,ui), £ k 'W +3 (<f> 2 ,u 2 ) < 2A. Define V(t) by 

m ~ 1 2TT1 " 02)|2 + ^ (ui " U2)l ^ 
2 r 1 

+ E / (2 fc + if - (V2m k h)\ 2 + lCVT)*(€*«x) - (^f (£V)1 2 ^ , 

where Vj, V^* are the V, induced by <f>j. We will prove that jgD < CP, which immediately 
gives the uniqueness. Let us consider i = 2 case only in V{t). Recall the system (jH]). 
By subtracting two systems from each other, we obtain the equations for (V*Vi(£ k <pi) — 



^ + 2(1 " 



n+1 
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V?V 2 (£ k <h), - V 2 V 2 *(t*u 2 )): 

dtWv^fa) - v;v 2 (t k <h)} - (2* + - v 2 v 2 \i k u 2 )} 

= (2k + l){v;v 2 V?(£ k U2) - V 2 *V 2 V 2 *(i k u 2 )} + {(Vi*)tVi(C fc 0i) - (V 2 *)tV 2 (C k h)} 
${W(£ fc «i) - V 2 Vm k u 2 )} + ^L^VxiVMefa) - V 2 *V 2 {i k <t>2)} 

' {V 2 V 2 *V 2 {e<l>2) - V x V2V 2 (£ k <h)} + 2{(V 1 ) t V 1 *(£ fc m) - (V 2 ) t V 2 *(Z k u 2 )} 



2k + 1 



(I) (II) 



We have to show that the L| norm of the right hand sides is bounded by T> 2 . We estimate 
(I) and (IT). Other two terms can be treated in the same way. First, we note that 



,2k+l j2fc+l 

K -9 2 



f\0 k {v l (i k 4>i) -v 2 {i k 4>2)}di' . 

Jo 



If simply applying Holder inequality, one gets 

,2k+l ,2k+l 

ll^p- - %rlUr - II^C* *i) - ^(eV 2 )IL| • (70) 
? 2 ? 2 

Apply Holder inequality once more: 

+ \\Vsv?(yi(Z k <h.)-V2(Z k fa))\\i%} 

w ' 

(**) 

Note that it is not trivial to get the boundedness of (**) in terms of T>?, since it is not of the 
right form yet. Here is the estimate of (**): 

(**)<ii^{wce fc 0i)- v 2 *v 2 (eh)}\y+\\vi{v 2 *v 2 (e<i>2)-v{v 2 (eh)}\\ 



j 2 



The second term can be written as 



Hence, by using ([70]) . we deduce that 



,2k+l ^k+l J2k ^k 1 

Note that one cannot hope to get the bound of ^fc+i ~~ ps+i with the P-regularity. Of course, 
it is bounded by A, but for the purpose of uniqueness, A-bound is not useful for the difference 
terms. The idea is to rearrange (I) and (II) so that P2 can be extracted from each term. 
We write the factors first and then L| factor at the last. 2? 2 can also come from L|° 
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factor thanks to the above estimate (|7ip . We start with (I): 

k2k ^2k A,2k 

I 



-(*) 

The first term is written as ) ^sf ^ 2 ^ 2 ~|F^ ^ 2 ' an d thus its L| norm is bounded by 

A and 2? 2. The second term can be treated as follows: 

~2fc + 10i £ £ fc 2 j £ fc+1 

Thus (*) is controlled by A and D2. Next we rearrange (//): 

2A: 0i C fc 01 02 ^ 

+ e e +l {e ^ [ 7T ]_e — e i — } 

2 



For i derivative difference terms, use the equation to convert into u terms and apply the same 
argument. For instance, the second term can be rewritten as 

e fe 0i e fe 1 01 02 j ' 

It is easy to deduce that (**) bounded by A and T>2 . Other terms can be similarly estimated. 
This finishes the proof of the uniqueness and Theorem 12.11 



7 Duality argument 

Here we would like to prove the existence for the linear problem ()60p . This is a consequence 
of the following proposition. 

Proposition 7.1. For f and g in L 1 (0,T;L 2 ) ; there exists a unique solution (F,G) to the 
linear system 

' d t F - V*G = f 

< d t G + VF = g (72) 
G{£ = 1) = 0, F{t = 0) = G(t = 0) =0 

on (0, T) which satisfies 

\\(F,G)\\ c{[0 , T] . L 2 ) <C\\(f,g)\\ LlL 2. (73) 
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Moreover, if (f,g) G L 1 (0, T;X fcj x Y ^) for some integer j < \k~\ and for < 2i < j — 1, 
we have {V*) 2i g = at £ = 1 and for 1 < 2i + 1 < j - 1, we have (V) 2i+1 f = at I = 1, 
then 

\\( F ,G)\\c(lo,T\;X'''ixY k >i) - C \\(f>9)\\mx k -3xY k 'i) ( 74 ) 
for some constant C which only depends on ||£ fc </>||xfc,rfcl+3 • 

The proof is based on a duality argument. Let A denote the set 

A = { rj) G C°°((0, oo) x (0, 1]), such that, = 1) = 0, (<f>, ^)t=T = 0}. 

Hence, (F, G) solves (|72p on a time interval (0, T) if and only if for each test function (</>, tp) G 
.4, we have 



-Fdt<f> - GVcp = I I f<j> 

(75) 



T r rT 



-Gdtip + FV*ij> = I gtjj. 
o J Jo J 

We denote = {^fc+VF) defined on the core 

\dt(i) ^LlLl (F,G)eL 2 (V(V),V(V*))}- 



Hence, C can be extended in a unique way to a closed operator. Moreover, A C £>(£*), the 
dual of £ and £*(*) = (~%£]Qf). Hence > flZSJ holds for each (<£,-0) G X if and only if for 
each (cj), ifi) G A, we have 

We take (J) G ,4 and denote £*($) = (*). The energy estimate written for C* yields 
that 



sup u\\h + u\\h<c f + wi 

0<i<T JO 



T 

11^112 , IUT.II 

L 2 ' 



Hence, the operator C* defines a bijection between A and C*(A). Let Sq be the inverse. 
Hence 

S : £*{A) -» A 



(?) ^ (I) 



and we have 



ll"5oQ)||c([0,T];L2) < C||^|| L i L2 . 



We extend this operator by density to C*(A) and to L X L 2 by Hahn-Banach. We denote 
this extension by S. Hence 

S: L X L 2 C([0,T];L 2 ) 
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Now, we want to solve ([72]) . namely = (^) with (g) i=0 = 0. This is, of course, 

equivalent to the fact that (fT6l) holds for each (4>, i/j) G A. 
Hence it is enough that for all (*) G L 1 L 2 , we have 



r 



(77) 



Therefore, it is enough to take (^) = 5* (/ g ) where S* is the dual of S which satisfies 

S* : M(0,T;L 2 ) -» L°°(0,T; L 2 ). 

In particular it maps L X L 2 into L°°L 2 . Hence (|73p holds with C([0,T];L 2 ) replaced by 
L°°([0,T];L 2 ). At this stage we do not know whether (^) is continuous. This will actually 
follow from the regularity. 

The uniqueness of (g) follows from the fact that if / = g = in (1721) , then (q) is the unique 
solution to (U2J| in L°°([0, T];L 2 ). To prove this consider a solution (g) G L°°([0, T];L 2 ) to 
(|72p with / = g = 0. We will also use the duality argument. Indeed, arguing as above and 
changing the roles of C and £*, we can prove the existence of a solution (t) G L°°([0, T]; L 2 ) 
to the dual problem 

-d t <f> + T/> = $ 

^ = i) = o, 

0(t = T) = ij}(t = 0) = 0. 
for each (*) G L l L 2 . 

Then, for each (*) G L^L 2 , we consider (^) a solution to (|78p . Hence, we can write (|76p 
with the solution (?). This yields 



r 



0. 



(79) 



which gives rise to the fact that F = G = 0. 

To prove (|74p. we can argue by induction of j. We start with the the case j = 1 and we 
first argue formally. Applying V and V* to (|72p. we get 



(80) 



Notice that the boundary condition VF(£ 
Hence, we deduce formally that 



d t V*G + V*VF = V*g + V t *G 
d t VF - VV*G = Vf + V t F 
VF{i = 1) = 0, 
[ VF(t = 0) = V*G(t = 0) = 0. 

1) = 0, comes from the fact that g = at £ = 1 



||(KF,V*G)|| 1 ,»< M . 12) < C\\(V'g + V t 'G, Vf + V,F)]\ L , V , 



To make this rigorous, we define first (^°), the solution of (|80p with the right hand side 



replaced by Cvf)- Hence, 



[yj!) ■ Then, we define for each integer i, (^) , the 
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solution of (j80j) with the right hand side replaced by Cvjr-iz* ^) = ). Hence 



\\(Y ,Z )\\ Laa{[0iTl . L2) < C\\(V*g,Vf)\\ L i L 2 
and using the fact that VtV^ 1 is bounded from L 2 to L 2 , we get that 

<C(CT)' || (F* 5 ,y/)|| LO o L 2 



Hence, denoting 



we see that ( z ) solves 



i=0 



<%y + y*z = v*g + y/y*-^ 
dtZ-yy = vf + v t v- 1 z 
z{i = i) = o, 
z(t = o) = y(t = o) = o. 



(81) 



The operator V can be thought of as an operator from T>(V) — ► L|. It has a transpose 
'V which goes from L| — > P(y)' where P(y)' is the dual space of T> (V). It satisfies : for 
/ G V{V) and u G L 2 , 

(V(f),u) L 2 iL 2 = (/,*y(«))©(y),27(V)'. 

Moreover, if we identify L 2 to a subspace of V(V)', then 'V extends y* to L 2 . 

The same argument shows that we can also extend V to L 2 by considering t V*. Hence, 
one can also interpret the system (|72p as equalities in V' and V*' and the same remark holds 
for dSH). 

Using that VldtiV^Z)] = d t Z - VtV~ x Z, we deduce that V[d t {V- l Z) - VY - Vf]. 
Since, the kernel of V is {0}, we deduce that 

d t (v- 1 z)-Y = f. 

We also have V^dtiV^Y)] = d t Y - VfV^Y. Hence, the first equation of flBU) can be 
written 

V*[d t (V*- 1 Y) + Z - g] =0. 

Hence since the kernel of V* is {0} when we add the vanishing of the boundary condition at 
£ = 1, we deduce that 

d t (V*- 1 Y) + Z-g = Q. 

By uniqueness for (|72p. we deduce that F = V~ X Z and G = V*~ l Y . Hence, we get that 
(F,G) G L°°(0,T; X k > 1 x Y^ 1 ). In particular, this also shows that {d t F,d t G G L l (0,T;L 2 ). 

Hence, integrating in time, we deduce that (F, G) G C([0, T); L 2 ). In particular this shows 
that (i73j) holds if have more regularity on (/, g). By a density argument, this shows that (f73|) 
holds even if we only know that (f,g) G L^l? . 

Arguing by induction on j, we prove that if (f,g) G L 1 (0,T;X k ^ x y^'- 7 ), then (F,G) G 



43 



C([0,T];X k 'i x Y k >i). 
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